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Among several analytic approximations for the growth of density fluctuations in the expanding 
Universe, Zel'dovich approximation in Lagrangian coordinate scheme is known to be unusually accu- 
rate even in mildly non-linear regime. This approximation is very similar to the Pade approximation 
in appearance. We first establish, however, that these two are actually different and independent 
approximations with each other by using a model of spheroidal mass collapse. Then we propose 
Pade-prescribed Zel'dovich-type approximations and demonstrate, within this model, that they are 
much accurate than any other known nonlinear approximations. 



When we analyze the growth of density fluctua- 
tions in the expanding Universe by analytical methods, 
the Zel'dovich-type approximations (ZTA hereafter) are 
known to be unusually accurate even in mildly non-linear 
regime for unknown reason |-|]. These Zel'dovich-type 
approximations are grounded on the Lagrangian coordi- 
nate scheme and are one-dimensional-exact; they become 
exact in the plain parallel mass distributions. The valid- 
ity of ZTA has been argued recently based on these phys- 
ical properties On the other hand, the appearance 
of these ZTA are very similar to the rational expansion 
method called Pade approximations. Though they have 
been widely used in the literature, the validity of these 
approximations has not yet been established as well. 

We would like first to compare these Zel'dovich-type 
and Pade approximation methods. It is almost impos- 
sible and is not pragmatic to argue in general analytic 
form. Therefore in this letter, we restrict our consid- 
eration to a model of spheroidal collapse which we can 
solve semi-analytically. If the above two types of approx- 
imations are actually the same, then they would give 
the same result for this restricted model. We demon- 
strate, in the first part of this letter, that this is not the 
case and conclude that they are independent approxima- 
tion schemes. Then this fact suggests a possibility to go 
beyond ZTA by Pade prescription on ZTA. We demon- 
strate, in the second part of this letter, that this Pade 
prescription dramatically improves ZTA. 

In the gravitational instability theory, the non- 
relativistic matter with zero pressure in an Einstein-de 



Sitter (EdS) universe is described by the following set of 
equations (see Rcf. fl0|]), 

5 + V ■[{\ + 5)v] = Q, (1) 
V + 2Hv + {v ■ \7)v + = 0, (2) 

= ^hH, (3) 

where x, v{x,t), ^{x,t) are respectively position, pecu- 
liar velocity, peculiar potential in comoving coordinate. 
The scale factor a varies as a oc i^/'^ and the Hubble pa- 
rameter is _ff = a/a = 2/(3i). Although we consider only 
EdS universe for simplicity in this letter, it is straight- 
forward to generalize the analyses in general Friedman- 
Lemitre universes. 

In the Eulerian coordinate scheme, the linear solution 
of them has a simple form Si^{x,t) = 6in{x)a{t) / ain, ne- 
glecting the decaying mode. Considering this to be a 
small parameter, we can naturally expand the full so- 
lution in powers of (5 = (5l + (J*-^-* + f^^'^' + • • •, 

$ = $L -f $(2) -I- $(3) _| , V = VL + t)'^' + v'-^'^ H , 

where (5*^"^ v'-"^ are assumed to be of order 

(e.g., see ||||). 

In the Zel'dovich-type approximations 0J|,^,|l^-^, 
we work in the Lagrangian coordinate scheme in which 
the location of a mass element x of the fluid is expressed 
by the initial location q and the time dependent displace- 
ment vector as a; (q, i) ^ q-\-'^{q,t). Then the density 
contrast S[x{q,t),t] = dct[dxi/dqj]~^ — 1 is determined 
by solving the equation of motion 
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where d/dt is the Lagrangian time derivative, J^^ 



dxi/dqj 



AeiJij. This Eq. 



is ob- 



tained from Eqs. (|^) and (^) , and Eq. (||) corresponds to 
the usual Eulerian vorticity-free condition. These nonhn- 
ear equations for * can be solved by the method of iter- 
ation considering d^i/dqj as small expansion parameter: 

•. In EdS universe, the time 



(3) 



dependence of each terms is separated from its spatial 
dependence: = (2/(3a2i72))n^(«) (-^^ rp^^^ 

order solution ijj'f''^ — —di(j)i,{q) is the original Zel'dovich 
approximation (ZA). 

Now we introduce a model of collapsing homogeneous 
ellipsoid. We parameterize the density perturbation 
5{x,t) as 



5(x,t)^5,{t) 6 1 



ai{t) ai{t) 



(6) 



where are the half-length of the principal axes of 
the ellipsoid and 6 is the step function. The solu- 
tion of the Poisson Eq. (||) inside this homogeneous el- 
lipsoid is known (see, Ref. |^6|) and it becomes $ = 

|7?^(5e X]i=i I ^^"i ttie equations of motion for 

three Ui are given by |p|,p7| 



(7) 



where 



The density contrast of the ellipsoid is obtained by ob- 
serving that 1 -I- (5 is inversely proportional to aiQ;2a3. 
These equations are solved by numerical integration. The 
solutions are 'semi-analytic' in this sense. In the follow- 
ing, we only consider spheroidal case, a\ — ai for sim- 
plicity. 

For the system of spheroidal perturbations, we apply 
ZTA, resulting in |9) 
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(23 + 23/iin - 39/if, 
- 25/i?„ + Uht 

- 46/iin - 39/if„ 
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30/ifj, (13) 



(14) 



where we have changed the parametrization, Ai = A2 = 
|(1 + h), A3 = |(1 - 2h) ||l|], and /ii„ = /i(tin), with 
iin being the initial time for numerical integration. The 
density contrast in these ZTA is given by 



S = 



919293 



{qi+<I'i){q2+1'2)iq3+1'3) 



-1, 



(15) 



where = ^^^'^ corresponds to the original ZA, = 
if-^i) -t-iZ^I^^ corresponds to post-Zel'dovich approximation 



(PZA), = tf/^^^ -hiZ'^^'' -htf'^"'' corresponds to post-post- 
Zel'dovich approximation (PPZA). 

In contrast to the above Lagrangian perturbation 
methods, the surface of the spheroid cannot be explic- 
itly expressed in Eulerian perturbation methods. We 
simply transform the expression already obtained in La- 
grangian perturbation scheme to that in Eulerian pertur- 
bation scheme. This is based on the fact that the small 
expansion parameters Wij in Lagrangian scheme and S in 



>(2) 



>(3) 



Eulerian scheme are the same order and thus ll'. 
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In our case, oc a", so Eulerian perturbative series 

can be simply obtained by expanding equation (^5|) in 
terms of expansion factor a: 



S = ±a+ { 1 ht 
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Now we introduce Pade approximation associated 
with the Perturbative expansions in Eulerian coordinate 
scheme. Pade approximation of type-(M, A'') for a given 
unknown function f{x) is expressed as the ratio of two 
polynomials. 



M 



fpadc{M,N) = I 



\k=0 




(17) 



The constant parameters ak and bk are determined so 
that they maximally yield the known Taylor expansion 
of the function f{x) Q up to (M + N)-th order. The 
density contrast in the spheroidal model up to the third- 
order perturbation is given by Eq. (|l^) . The correspond- 
ing Pade approximation of type-(l,2) is given by 
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We observe that the ZTA density contrast of Eq. ( p^ ) 
and Pade density contrast of Eq. (|l^) are definitely dif- 
ferent with each other. Actually the plots of these ap- 
proximations in Fig. 1 for spherical perturbations clearly 
demonstrate the difference. Some of lines in Fig. 1 have 
been previously appeared |^,|| . 
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we conclude our work. First we focused on th: 



th it ZTA am.d Pade approximations are actually 
scheme orsnot. By using spheroidal pertilrbal 
L, we found tj^at these two types of approxiijiaj^i 
different indepeh^ent scheme with each other (^ei 
appearanc^ Then, based on this fact, we' 
the Pade prescription for the ZTA. We found 




el'dovich, Astron. Astrophys. 5, 
^el'dovich, Astrophysics 6, 164 ( 
shi, V. Sahni and A. A. Starobii 
517 (1994). 
and P. Coles, Phys. Rep. 262, 
and S. Shandarin, Mon. Not. 
(1996). 

A. L. Melott and A. G. W( 
. 288, 349 (1994). 
Melott, T. Buchert and A. G. We 

294, 345 (1995). 
Bbuchet, S. Colombi, E. Hivon and 
Astrophys. 296, 575 (1995). 
lisato, T. Matsubara and M. M 



■o-ph/9707296 



Peebles, The Large-Scale Struc 
ehon Univ. Press, Princeton, 1980) 
': y, Astron^s. J. 279, 499 (1984) 
Cjlpff,^^ Grinstein, S.-J. Rey a 
AstrcJpJ|iys. J. 31^6 (1986). 

and^:^hlers, Mon. Not. R. A 



B4cl.ert 



1Q93). 

lert Mon. Not. R. Astron. Soc. 2 
ate Ian, Mon. Not. R. Astron. Soc. 2 
)gg. Foundation of Potential The 
395). 

. White and J. Silk, Astrophys. . 
efficients Ai automatically satisf; 
2, which follows from Poisson E 
Render and S. A. Orszag, Advanc 
for Scientists and Engineers (M 
Compaiy, New York, 1978). 



00 

o 



d 



enj} , enj} xojddB 

■0-= 9] 9/ M 



